Metal-insulator transition and nature of the gap in NiS2-xSe a; 
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The origin of the gap in N1S2 as well as the pressure and doping induced metal- insulator transition 
in the NiS2- a; Se :!; solid solutions are investigated both theoretically by the first-principles bandstruc- 
tures combined with the dynamical mean-field approximation for the electronic correlations and 
experimentally by means of infrared and X-ray absorption spectroscopy. The bonding-anti-bonding 
splitting in the S-S (Se-Se) dimer is identified as the main parameter controlling the size of the 
charge gap. The implications for the metal-insulator transition driven by pressure and Se doping 
are discussed. 

PACS numbers: 71.30.+h, 62.50.-p, 78.30.-j 



I. INTRODUCTION 

The metal-insulator transition (MIT) due to electronic 
correlations has been subject of intense research for sev- 
eral decades. The NiS2- a; Se 3 ; series provided an impor- 
tant model system exhibiting a MIT controlled by vary- 
ing the Se content x, temperature T or pressure P— ~— . 
Particularly interesting is the similarity of the x — 
and P — T phase diagrams^ to that of the Hubbard model 
in the infinite dimension limit 6 , consisting in the presence 
of the phase transition between paramagnetic metal and 
paramagnetic insulator at intermediate temperatures, ex- 
istence of the high-temperature crossover regime, and the 
possibility of driving a metal insulating, in a certain range 
of x, by increasing the temperature. Consequently, the 
MIT has been commonly attributed to broadening of the 
Ni-d band 3 - - — . Despite a large volume of available 
experimental data, the microscopic origin of the MIT 
in NiS2- a; Se a ; is poorly understood and a satisfactory 
material-specific theory is missing. 

NiA 2 (X=S,Se) can be viewed as NiO with the O 
atom replaced by an X2 dimer. Strong hybridization 
between the X — p orbitals pointing along the X — X 
dimer, p a orbitals, leads to a formation of split bond- 
ing and anti-bonding bands (the latter ones are referred 
to as p*), which accommodate two holes leading to 
an X^ 2 valence state. Numerous photoemission (PES) 
studies 7 - ' 10 ' 11 complemented by cluster calculations^ re- 
vealed a similarity between the NLY2 and NiO valence 
band (VB) spectra, supporting this analogy. The band- 
structure calculations in the local density approximation 
(LDA)^! 2 *! 3 . rendered NiS 2 a metal with a partially filled 
e g band and empty p* orbitals. By analogy to NiO it was 
speculated that local d — d correlations open a gap be- 
tween the S-p band and the upper Hubbard band of Ni-d 



e g character—, leading to the classification of MS2 as a 
charge-transfer (CT) insulator in the Zaanen-Sawatzky- 
Allen (ZSA) scheme 1 ^, although in the original ZSA pa- 
per NiS 2 was classified as 'p-type metal were it not for 
the gap between the S p% and the rest of the p band'. 

We have combined theoretical calculations with new 
experimental data in order to address the following ques- 
tions: i) Why is NiS 2 an insulator and NiSe2 a metal?, ii) 
How do N1X2 respond to external pressure?, iii) Are the 
effects of pressure and varying Se content equivalent?, 
iv) What is the mechanism of the MIT in NiS2- £C Se a: ? 
The electronic structure of N1A2 is studied numerically, 
using a combination of the ah initio bandstructure and 
the approximation of the dynamical mean-field theory 
(LDA+DMFT)&±£. The infrared (IR) reflectivity is mea- 
sured on samples with various Se concentrations under 
applied pressure to investigate the evolution of the op- 
tical gap and the spectral weight transfer. The x-ray 
absorption spectroscopy (XAS) at S K-edge is used to 
elucidate the orbital character of features in the single- 
particle spectrum of NiS 2 . 



II. COMPUTATIONAL METHOD 

Our computations proceed as follows. First, we per- 
form paramagnetic LDA calculations 1 - 6 - and minimize 
the total energy to find the equilibrium S(Se) positions. 
Next, we represent the one-particle Hilbert space of 
the hybridized Ni-d and S(Se)-p bands in Wannier ba- 
sis (44 orbitals per unit cell) and calculate the parame- 
ters of the on-site d — d interaction with the constraint 
LDA approac h 17 ' 18 . We found only moderately different 
C^Nis 2 = 5 eV and UuiSe 2 = 4.7 eV and used J of 1 eV 
throughout the study. Finally, we construct a multi-band 
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FIG. 1: (color online) The crystal structure of N1S2: Ni atoms 
red (bright), S atoms blue (dark). The S-S dimer is high- 
lighted in the center of the figure. 

Hubbard Hamiltonian 

H = alH k a k - e dc (n d )N d + £ hf U dd hf , (1) 

k ieNi 

and compute the single-particle spectra using the DMFT 
approximation. Here, is the LDA Hamiltonian, 
44x44 matrices (four N1X2 units) in the basis of Ni-d and 
X-p Wannier orbitals, on a uniform fc-mesh in the first 
Brillouin zone^. The second term is the double-counting 
correction amounting to a constant shift applied to Ni-ii 
site energies^. The last term is the two-particle inter- 
action at the Ni sites in the density-density approxima- 
tion. The DMFT equations are solved on the Matsubara 
contour. The continuous time quantum Monte-Carlo im- 
purity solver— is employed to solve the auxiliary multi- 
orbital impurity problem. Single-particle spectral densi- 
ties are obtained by analytic continuation of Monte-Carlo 
data using the maximum entropy method 2 ^. 

III. EXPERIMENTAL DETAILS 

Single crystals, grown by vapor transport technique, 
were slightly polished in order to obtain a clean, mirror- 
like surface. Near-normal incidence reflectivity was mea- 
sured with a cassegrainian-based IR microscope coupled 
to a Michelson interferometer. The NiS2-xSe K series 
(x=0, 0.3, 0.5, 0.7) was studied at ambient conditions, 

A clamp diamond anvil cell (Diacell cryoDAC-Mega) 
equipped with type Ha diamonds was employed to per- 
form pressure dependent studies. A small piece of MS2 
(NiSi.7Seo.3) was cut and placed together with Csl (as 
hydrostatic medium) inside a 400 /im hole drilled in a 



90 /im thick CuBe gasket. The pressure was evaluated 
in situ, by the standard ruby-fluorescence techniquei2i 
Measurements were carried out partly at the University 
of Augsburg and partly at the IR1 beamline in ANKA, 
so to exploit the high brilliance of synchrotron radiation 
in the far-IR range. We recall that due to the small 
size of the samples inside the diamond anvil cell mea- 
surements as a function of pressure were performed for 
uj >250 cm -1 . The measurement procedure, as well as 
the data processing to extract the o ptic al conductivity 
oi(uj) are carefully described in Refs. 124051 

The sulfur Is X-ray absorption spectra (XAS) spectra 
which provides information about vacant S 3p-states were 
measured at the soft X-ray microcharacterization beam- 
line (SXRMB) at the Canadian Light Source. The spec- 
tra were measured in total electron yield (TEY) mode. 
Additional measurements were taken in total fluorescence 
(TFY) mode using a channel-plate detector, and partial 
fluorescence (PFY) mode using a silicon drift detector 
with a resolution of about 100 eV. Both the PFY and 
TFY measurements showed the same features as the TEY 
measurements, indicating that the surface of the samples 
was clean. The TEY measurements were used for the 
analysis because they had a better signal-to-noise ratio 
than the TFY or PFY measurements, and did not suf- 
fer from the self- absorption effects present in the fluores- 
cence measurements. An Si(lll) crystal was used in the 
monochromator, and the experimental resolution (E/E) 
was about 10000. 



IV. RESULTS AND DISCUSSION 

We start the exposition of our results with a brief re- 
view of the non-interacting LDA bandstructures, which 
agree with those published previousl y 8 ' 12 ! 13 . The Ni-d or- 
bitals form narrow ti g and broader, Ni-X anti-bonding, 
e a g bands due to octahedral coordination of Ni with the 
chalcogen X. The deviations from an exact octahedral 
symmetry at the Ni sites lead to further splitting of 
the tig bands into an doublet and an a\ g singlet. 
This splitting plays no discernible role in our investiga- 
tion. The characteristic feature of the chalcogen bands is 
a large bonding-anti-bonding splitting of the p-orbitals 
pointing along the X-X dimer (Fig. [IJ . Due to the S-S-Ni 
angle being close to the right angle («104 deg in MS2) 
the dimer states hybridize only weakly with the Ni e a g 
states. 

The Ni-X and X-X dimer distances as function of the 
lattice constant were obtained by minimization of the 
LDA total energies (see Fig. Q}. At ambient pressure 
the calculated MS2 bond lengths agree well with exper- 
imental observations. Compression leads to shortening 
of the Ni-S distance while the S-S dimer behaves as a 
rigid object. The situation is less clear in NiSe 2 . While 
the calculated ambient pressure Ni-Se bond length re- 
produces the experimental value, the length of the Se-Se 
dimer is substantially overestimated. Under pressure our 
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FIG. 2: (color online) The calculated pressure dependence of 
the equilibrium Ni-X and X-X distances. The isolated points 
mark the experimental values of Refs. [ll| (star) and[27] (cross). 



calculations yield an initial increase of the Se-Se indicat- 
ing that the Ni-Se bond becomes stronger relative to the 
Se-Se bond and suggesting that Se-Se dimer cannot be 
viewed as a molecular unit. We will discuss later on how 
the behavior of the bond lengths affects the electronic 
properties. 



A. Single-particle spectra 

1. MS 2 

In this section we present orbitally resolved spec- 
tral densities obtained at ambient pressure with the 
LDA+DMFT method. In Fig. [3] we show the NiS 2 spec- 
tra, which exhibit a small charge gap in agreement with 
experimental observations. The valence band (VB) re- 
sembles that of Ni02& with the characteristic distribu- 
tion of the Ni-d spectral weight between d 7 (around -6 
eV) and d 8 L excitations (at low binding energy), the re- 
sult of a strong p — e g hybridization. Overall, a good 
correspondence between positions of the main features 
in the experimental x-ray photoemission spectra (XPS) 
and their theoretical counterparts can be stated. Partic- 
ularly important for the further discussion is the high- 
energy shoulder at -7 eV, which is dominated by the S-p 
density in agreement with the observation of Ref. [ijj. 
The absence of an obvious gap in the experimental spec- 
tra shown in Fig. [3] is due to the experimental resolution. 

The calculated conduction band (CB) spectrum con- 
sists of two peaks: (a) p* band with a large resonant e g 
contribution at 1 eV, and (b) the upper Hubbard band 
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FIG. 3: (color online) The calculated ((7=5 eV, T=580 K) 
orbitally resolved spectra of N1S2 [Ni-d:eg (black) and ti g =e ll g 
(red) + ai g (red - dashed); S-p: total (blue, thick) and p a 
(blue-shaded)] compared to the experimental XPS+BIS (gray, 
shaded)^. The a, b, and c are denned in the text. The inset 
shows the conduction band e a g spectra for (7=5, 6, 7, and 
8 eV. 



of pure e g character at 3.5 eV. The identification of b 
with the upper Hubbard band is based on the observa- 
tion of its shift upon increasing the interaction strength 
(see the inset of Fig. [3]) and the absence of p charac- 
ter. The experimental CB spectrum, obtained with the 
Bremsstrahlung isochromat spectroscopy (BIS) of Ref. 
[Til , exhibits three distinct features: a sharp low-energy 
peak at 1 eV, a high-energy peak at 3.5 eV and a broad 
band (c) starting around 5.5 eV. The 5.5 eV band is 
not spanned by the Hamiltonian (JXJ) and thus the corre- 
sponding contribution is missing in Fig. [3] The orbital 
character of the other two peaks is of key importance for 
understanding the nature of the charge gap. The ques- 
tion can be formulated simply as: 'is the upper Hubbard 
band located below or above the p* band' ? The authors 
of Ref. 0j] interpreted the 1 eV as Ni-d upper Hubbard 
band and the 3.5 eV peak as S p* band, which is a clear 
contradiction to our result. 

In order to resolve this issue we have measured the 
XAS at the K-edge of sulfur, which can be directly com- 
pared to the calculated S-p spectral density as shown in 
Fig. |U The theoretical high energy spectrum was ob- 
tained by projecting out the S-p contribution from the 
LDA bands not included in the Hamiltonian ([T]). The 
identification of peak a with p* is corroborated by its po- 
sition with respect to the X-ray emission (XES) spectrum 
as well as by its distance from the high energy (5 eV) XAS 
band. Although the calculated onset of the high energy 
band is misplaced by about 1 eV with respect to the ex- 
periment, the overall agreement is very good. We point 
out that the experimental XAS a and c peaks match their 
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FIG. 4: The calculated S-p spectral density obtained from 
LDA+DMFT (full line), augmented with the S-p contribution 
from the high energy LDA bands (dashed line), compared 
to experimental XES spectral of Ref. |28| (circles) and XAS 
(diamonds) spectra of this work. 
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FIG. 5: (color online) The calculated ([7=4.5 eV, T=580 K) 
orbital resolved spectra of NiSe2 compared to the experimen- 
tal XPS+BISAi spectra, (the same notation as in Fig. [3| The 
position of the calculated high band is marked with the dotted 
line. The left inset shows the e a g densities for various interac- 
tion strengths U (eV). In the right inset the effect of reducing 
the Se-Se distance (A) on the spectra (C/ =4.5 eV) is shown. 



2. NiSe 2 

The LDA+DMFT calculations yield NiSe 2 metallic. 
Its low-energy spectrum is rather sensitive to changes 
of U, as shown in the left inset of Fig. [5] For U <5 eV, 
a temperature-dependent peak appears at the chemical 
potential. Nevertheless, the metallic state is robust and 
the spectrum remains gapless up to the largest studied 
value of U=8 eV (not shown). 

The VB spectra NiSe2 (Fig. [SJ overall resemble the 
N1S2 ones with some quantitative differences, such as the 
shift of the high-energy shoulder to -6.5 eV, also visible 
in the experimental data. As in N1S2 we find the p* 
band at the bottom of the CB manifold. Comparison 
of the calculated CB spectrum with the BIS data is less 
satisfactory than in the case of N1S2. In particular, the 
experimental 2 eV peak is not found in the calculation. 

Bearing in mind that the Se-Se dimer length of 2.52 A 
is overestimated, we have performed a reference calcula- 
tion with a reduced dimer length of 2.37 A (just below 
the experimental values). As shown in the right inset of 
Figf5J an increase of the dimer bonding - anti-bonding 
splitting leads to a strong suppression of the spectral den- 
sity at the lowest energies (see the right inset of Fig. [5]) , 
but is not sufficient to open a charge gap. Since the 
bonding-anti-bonding splitting between the VB bottom 
and the p* peak is essentially a molecular property of the 
Se-Se dimer, it does not depend on the value of U. 



B. Origin of the charge gap 

To analyze the behavior of the p* band we have calcu- 
lated the bare p bandstructures, obtained by decoupling 
the Ni-d states. This can be done conveniently by re- 
moving the Ni-centered orbitals from the LDA basis set 
[Fig. 23c), (d)], or rigorously using the p — p block of the 
Hamiltonian in Wannier basis [Fig. E)Ja),(b)] with quali- 
tatively similar results. The p-bands form a broad man- 
ifold bounded by dimer bonding bands at the bottom 
and a distinct anti-bonding p* band complex at the top, 
separated by a (pseudo)gap. The key difference between 
NiS2 and NiSe 2 is the size of this gap separating the p* 
band from the rest of the p-manifold controlled by the 
length of the dimer. The gap in MS2 is large enough to 
sustain the coupling to the correlated Ni-d orbitals, while 
the NiSe2 gap is too small and the coupling to d orbitals 
closes it completely. 

Next, we consider the effect of applied pressure. The 
LDA+DMFT calculations for pressures up to 13 GPa 
show the closing of the charge gap (Fig. [7]). This is a 
result of the combination of the rigidity of the S-S dimer 
and an overall band broadening. As a consequence of the 
former the bonding - anti-bonding splitting is insensitive 
to pressure. On the other hand, the widths of the in- 
dividual p-bands, controlled by the reduced inter-dimer 
distances, increase with pressure thus reducing the bare 
p gap. Together with increasing p — d hybridization this 
leads to the metal-insulator transition. Noteworthy the 
LDA results showed an elongation of the Se-Se dimer at 
moderate pressures, suggesting that the pressure-induced 
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FIG. 6: (color online) The bare S(Se)-p spectral densities 
and bandstructures. The spectral densities in the left panel 
[(a),(b)] are obtained from a Wannier function Hamiltonian, 
while the bandstructures in the middle panel [(c), (d)] are 
based on a LAPW calculation with removed Ni-like basis func- 
tions. Dashed bands for NiSe2 correspond to a reduced Se-Se 
distance. 



closing of the charge gap proceeds faster in compounds 
with higher Se content. 



C. IR reflectivity 

To verify the theoretical conclusions and to obtain 
further information about the electronic properties of 
NiS2- x Se a we have measured the infrared properties of 
the NiS2- a; Se 2: series (x=0; 0.3; 0.5; 0.7) at ambient con- 
ditions, shown in Fig|5k). For N1S2, R(w) slightly in- 
creases with increasing frequency, showing a broad max- 
imum in the mid-infrared and two narrow phonon ab- 
sorption in the far-IR. As Se is chemically substituted 
in the compound, R(w) increases its absolute value at 
low frequency, partially shielding the phonon lines and 
eventually becoming metallic-like (i.e. decreasing with 
increasing u>) for x= 0.7. The corresponding optical con- 
ductivities (j 1 (co) clearly show how the transition from an 
insulator (N1S2) to a metal (NiSi.aSeo.?) takes place: A 
well-defined gap followed by an absorption band is found 
for N1S2 . As Se is introduced the absorption band moves 
towards lower energies, progressively closing the gap. 
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FIG. 7: (color online) The comparison of N1S2 spectra in 
the insulating state at ambient pressure (upper panel) and 
in the metallic state under compression (lower panel). The 
black line marks the e^, is shaded area (blue) is the total S p 
contribution. The i2 9 spectra which show very little change 
with pressure are not shown. 

A similar behavior is found as pressure is applied: as 
shown in FigJSfc) it is again the red-shift of the absorp- 
tion band that closes the gap. A closer inspection at the 
curve at 5 GPa for MS2 shows the coexistence of two 
distinct terms at low-frequency. As it was also shown in 
recent IR measurements^, our data can be fitted at low 
frequency with two distinct features: a Drude-like term 
and a Lorentz peak. 

Our aim is to follow and compare the gap-closure by 
chemical alloying and external pressure. We define the 
experimental optical gap as the frequency at half-height 
of the absorption band and report its values in Figj9] 
for MS2 and NiSi.7Seo.3 under pressure. This defini- 
tion of the gap leads to a similar pressure dependence 
as obtained by linearly extrapolating the steepest rise of 
u\ (oj)^ or by considering the maximum of the absorption 
peak obtained with Drude-Lorentz fitting^. We remark 
that strong phonon absorptions of the diamond prevent 
to have (fully) reliable data between 1700-2500 cm" 1 . 
To overcome this problem and unambiguously define the 
gap value, the data analysis was performed as follows: 
reflectivity data are fitted with a Drude-Lorentz model 
in order to extrapolate the missing parts of the spectra. 
The KK transformations are performed, followed by a 
simultaneous fitting of R and sigma. This process is car- 
ried out iteratively until the <j\(lu) is stable with respect 
to the extrapolation used in the KK transformations. 
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FIG. 8: (color online) (a) Room-temperature optical conduc- 
tivity spectra of NiS2- a: Se a : (x— 0; 0.3; 0.5; 0.7, as indicated 
in figure) at ambient pressure. In the inset the reflectivity 
R(w) is reported, (b) Optical conductivity of N1S2 (solid) 
and NiSi.7Seo.3 (dashed) for increasing pressure. 
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FIG. 9: (color online) (a) The gap vs pressure for NiS2 and 
NiSi.7Seo.3 (note the mutual horizontal offset) 3(c) SW trans- 
fer ASW for NiS 2 at 5 GPa and for NiSi. 3 Se .7 at ambient 
pressure. 



Comparing the conductivity spectra we find that 
the qualitatively similar effect of pressure on NiS 2 
and S substitution by Se follows the scaling relation 
P[GPa]«x/0.15 proposed in Ref. 0. The pressure depen- 
dence of the optical gaps for pure NiS 2 and NiSi.7Seo.3 is 
shown in Fig]^), with an horizontal off-set according to 
the scaling relation. Our experimental results show that 



in the Se-doped compound pressure is more efficient in 
closing the gap, in accord with the previous theoretical 
reasoning. 

Next we discuss the metallic phase. To quan- 
tify the 'metallicity' we use the spectral weight 
(SW) transfer relative to the ambient pressure 
NiS 2 data ASW=SW(x, P)-SW(at=0,P=0), where 
SW(w)= J Q a\(bj')dbj' . Recalling that at low frequency 
there is the coexistence of two terms we can infer 
that with applied pressure a large part of the spectral 
weight is transferred to the mid-IR region first, while 
the appearance of Drude-like peak occurs upon further 
pressure increase [see Fig. [EJb)]. Without making a 
quantitative connection we point out that a qualitatively 
similar behavior is observed in the DMFT results: As 
pressure is applied to NiS 2 , due to the rigidity of the 
S-S dimers the transition occurs mainly due to band 
broadening not corresponding at first instance to the 
formation of a Drude term. Instead, when the Se-Se 
distance is varied we expect a fast response of the low- 
energy electrodynamics (right inset of Fig. [5]). Following 
the scaling relation we compare NiS 2 at 5 GPa with 
NiSi.3Seo.7 at ambient pressure in Fig. (9b). Although a 
similar amount of SW is transferred below 4000 cm" 1 , 
NiSi.3Seo.7 exhibits a higher SW transfer at lower 
frequencies (<2000 cm -1 ) pointing to a more intense 
Drude-like term, also supporting our conclusion about 
the faster onset of metallicity with Se substitution. 



V. CONCLUSIONS 

Using the LDA+DMFT approach we have reproduced 
the paramagnetic phases of NiS 2 (insulator) and NiSe 2 
(metal) and provide the following answers to the ques- 
tions posed in the introduction: i) The presence of the 
charge gap in NiS 2 in contrast to NiSe 2 is a consequence 
of a larger gap in the bare p-bandstructure (p-gap) of 
NiS 2 resulting from stronger bonding within the S-S 
dimer. ii) The pressure reduces the p-gap by broadening 
the individual p-bands, while the Se substitution reduces 
the p-gap due to smaller bonding-anti-bonding splitting, 
iii) The MIT is controlled by varying the size of the p-gap. 
We find both computationally and experimentally that 
the pressure driven MIT is accelerated by the Se substi- 
tution. The conceptual picture of NiS 2 _ x Se x is provided 
by a periodic Anderson model with a variable gap in the 
conduction electron bath rather than the charge-transfer 
insulator model. 
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